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ABSTRACT 

We extend the approximate radiative transfer analysis of Hershkowitz, 
Linder, and Wagoner (1986) to a more general class of supernova model 
atmospheres, using a simple fit to the effective continuum opacity produced 
by lines (Wagoner, Perez, and Vasu 1991). At the low densities considered, 
the populations of the excited states of hydrogen are governed mainly by 
photoionization and recombination, and scattering dominates absorptive 
opacity. We match the asymptotic expressions for the spectral energy density 
Jy at the photosphere, whose location at each frequency is determined by 
a first-order calculation of the deviation of Ji, from the Planck function By. 
The emergent spectral luminosity then assumes the form Ly = 4'n''^rl('^By{Tp), 
where Tp{u) is the photospheric temperature, ( is the dilution factor, and is 
a fiducial radius [ultimately taken to be the photospheric radius Tplu)]. The 
atmosphere is characterized by an effective temperature Tg (oc L'^^'^r* ^^^) and 
hydrogen density uh = n^,{r^,/r)"' ; and less strongly by the heavy element 
abundance and velocity gradient. Our major result is the dependence of C on 
frequency z/ and the parameters Tp, Vp, and a. 

The resulting understanding of the dependence of the spectral luminosity on 
observable parameters which characterize the relevant physical conditions will 
be of particular use in assessing the reliability of the expanding photosphere 
method of distance determination. This is particularly important at cosmological 
distances, where no information about the progenitor star will be available. 
This technique can also be applied to other low-density photospheres. 
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1. Introduction 

The goal of this investigation is to gain a broader understanding of the dependence of 
the continuum spectral luminosity of Type II supernovae on the physical conditions near 
their photospheres as reflected by parameters that are in principle observable. We extend 
the semi-analytic approach of Hershkowitz, Linder, and Wagoner (1986b, hereafter HLW) 
to a more general class of Type II supernova model atmospheres, using a simple fit to the 
effective continuum opacity produced by lines (Wagoner, Perez, and Vasu 1991; hereafter 
WPV). Understanding the continuum spectrum allows us to calculate the dilution factor, 
a frequency dependent quantity that appears as a correction to the luminosity due to 
the non-blackbody, reduced flux of the continuum spectrum emitted by supernovae. This 
dilution factor is arguably the most critical element in the expanding photosphere method 
(EPM) of distance determination (Kirshner & Kwan 1974; Wagoner 1981; Eastman & 
Kirshner 1989; Schmidt, Kirshner, & Eastman 1992, hereafter SKE; Wagoner & Montes 
1993; Eastman, Schmidt, and Kirshner 1994; Schmidt et al. 1994). 

The extensions of this work from that of HLW include: 1) spherically-symmetric power 
law atmospheres; 2) inclusion the UV opacity of the heavy elements via a simple fit from 
WPV; 3) allowance for atmospheric regions of two different types: mostly ionized {ue = uh', 

case I), and mostly neutral (n,. <^ O.Sn^; case N). In addition, at the low temperatiues 
studied (5, 000 ^Tp ^ 20, OOOK) it is often found that the Lyman continuum is in radiative 
detailed balance. Then collisional coupling of the n = 1 and n = 2 levels of hydrogen 
becomes the most important channel determining the departure coefficient of the ground 
state. Thus we retain collisional coupling between these levels, although we have found that 
it is otherwise negligible at the low densities of these photospheres {uh <^ 10^^ cm~^). 

We would like to stress that although our results are approximate, they require no 
knowledge of the progenitor star and are expressed in terms of observable parameters which 
fully characterize the relevant physical conditions near the photosphere. The analytical 
aspect of our results also reflects the pedagogical goal of our investigation. 
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2. Radiative Transfer Problem 

The assumptions which define our radiative transfer problem are hsted below. 

1. The velocity gradient is small enough to yield a quasi-static photosphere. This is a 
good approximation for most supernovae after a few days. 

2. The supernova is spherically symmetric, with the total hydrogen density given by 
riH — n^{r^/r)'^ near the photosphere. Here and r* are fiducial values. 

3. We consider hydrogen dominated systems but include heavier elements via an 
effective continuum scattering opacity (from WPV). Hydrogen photoionization and inverse 
bremsstrahlung account for the absorptive opacity. 

4. We consider only low density photospheres, such that collisions may usually be 
neglected (except for the ground state of hydrogen) and for the optical and UV photons the 
atmosphere is scattering dominated (by electrons and mostly iron-peak lines). 

5. Radiative detailed balance holds for the hydrogen lines, as verified by Hershkowitz 
and Wagoner (1987). 

6. We shall adopt the value of 1/3 for the Eddington factor near (and below) the 
photosphere since in scattering dominated photospheres the continuum is formed at an 
optical depth large enough to make the radiation field essentially isotropic. 

7. Energy is solely transported by radiation. 

8. In keeping with our approximate analysis, the hydrogen Gaunt factors are set to 
unity. 

Assumptions 2,3 & 4 are more general than those made in HLW, who also only 
considered the case of complete ionization and did not examine the dilution factor in any 
detail. The total opacity and optical depth are 

roo 

Xv = Xsc + i^u, Xu{i^, r')dr' , (1) 

Jr 

where the scattering opacity Xsc — Xes + Xunes- The absorptive opacity is due to hydrogen 
photoionization and inverse bremsstrahlung: 

= E(^^ - n*e-^^/'^)a,M + ^e«//(^, T){1 - 6'^^/'^) . (2) 

i 

The LTE populations are given by n* — nl^i{T), where ^i{T) is the Saha-Boltzmann 
function for a pure hydrogen gas. We assume (and verify) that most of the free electrons 
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come from hydrogen, giving 

nH = ne + {l + di)nl , (3) 

where di is the ground-state departure coefficient. 

For our extended atmosphere, we use the Mihalas (1978) formulation of the transfer 
equations: 

d{UqyJu)/dT,, = U/iAnf , (4) 
dL^ (47rr)2 



[KyJu - nlB^) . (5) 



In these equations fi, is the Eddington factor, drg^, = q^dr^, is the 'spherical optical depth', 
and -Bjy = {2hi'^/c'^)[exp[hi'/kT) — 1]^^ is the Planck function. The sphericality factor q^y 
is given by ln(r^g,^) = /J'J(3/i. — 1) / {r' fi^)]dr' + Inr^. Since for scattering dominance the 
total optical depth at the photosphere rp(z/) > 1, we may use fu{T,y ^ Tp) = 1/3, giving 
qu = {j'^:/r)'^. Equation is then simplified to: 

dJy SLi^ , , 

(6) 



dTgy (47rr* 



|2 



The constraint of radiative equilibrium imposes the equivalent conditions 

/•oo 

/ L,ydiy = AnrlajiT^ = constant , (7) 
Jo 

POO 

/ (k.J, - <5,)c/z/ = . (8) 
Jo 

The effective temperature Tg and fiducial radius r^, are parameters characterizing our model 
photospheres. The second constraint is that the atomic level populations are maintained 
in statistical equilibrium (not LTE). Because of our assumption of relatively small velocity 
gradient broadening, we can take the radiative transitions between the bound states to 
be in detailed balance (Hershkowitz & Wagoner 1987). That is, there is a large optical 
depth in the lines at the continuum photosphere. Then the non-LTE departure coefficients 
di = {ni — nj) / rii are given by 



dj 



(■oo ^ roo 

An / {ai^J,/hu)du + ^ - ^ dfii, = An / {B, - J.) [1 - e"^''-"'^^] {aijhv)dv , 
Ju^ J Ju, 

(9) 

where niCij is the rate of transitions from level / — > j induced by electron collisions (Mihalas 
1978). [Here n denotes unbound states, ui are the threshold ionization frequencies, and L 
(usually taken to be 10)is the principal quantum number of the highest level included.] 
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The appropriate second order equation is obtained by substituting equation (||) into 
equation (13) , which yields 



[K^Jy K By) . 



(10) 



We rewrite the bound-free and free-free absorption coefficients (setting the Gaunt factors 
to unity) as 



(11) 
(12) 



The constraint (^ of statistical equilibrium can be incorporated directly into the radiative 
transfer equation (0), yielding for z/„ < z/ < z/„_i our working transfer equation 



3ne 



AffiT)+j:MT) 



l=n 



l=n 



Here M^^- ^ is the inverse of the matrix 



K13) 



which is constructed from the relative transition rates 



An / {ai^Jy/hij)diJ 



-1 



(14) 



(15) 
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3. Approximate Analysis 



We proceed by first matching asymptotic solutions to the transfer equation at the 
photosphere. This gives an expression for the emergent luminosity which depends upon the 
location of the photosphere. The location is then determined by solving equation ([131) to 
first order in | Jj, — B^\/B^. Then the relevant physical quantities at the photosphere can be 
related. 

For ^ Tpi'^) we recover LTE populations and a thermalized radiation field, giving 
Jy = By. Integrating equation (^) over frequency, we obtain the familiar result 

= 16^ + ^) " -J^^'" + C)AC = const) . (16) 

Here is the Rosseland mean spherical optical depth [dr^ = —qxRdr). In the other 
asymptotic limit -C Tp(i/), the spectral luminosity is no longer changing {Ly = L°), 
although scattering will still be important where Ty 1. Integrating equation (^) in this 
regime yields 

= 3L°(47rr*)"^(rg,, + Cy) , {cy = const.) . (17) 

If we now match the asymptotic expressions for Jy at Ty = Tp{v) we obtain the emergent 
luminosity 

LI = A7rrlC\By{Tp) , (18) 

where the dilution factor ( is given by ("^ = {4/3){Tqy + Cy)^^. Near the photosphere we take 
Tqu = {,Xy/XR)'^q > ~ C ~ 1, giving 



1 



Xu I „ \T, 



p 



p, 



The subscript p indicates that the quantity is evaluated at the (frequency dependent) 
photosphere. 

In order to proceed, we must approximate the dominant opacities in various portions 
of the spectrum. Employing the results of WPV, we adopt the form 

1 + no$iaiK/(Tes if z/ > ULyman and 9 <9eg , 

1 + {no^iY^'^aii^/aes if z/ > uiyman and 6 > 6eq , 

^=1 1 + (250^^16 + m^){iy/l^Lymany/^ if l^Lyman > > l^inter , (20) 
1 + (250^16 + 10^4) (^/^.^^^^)10 if ^.^^^^ > ^ > Upaschen , 

1 otherwise , 
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where Xes = neCJes ,0 = 5040/T(K) and I'mter = c/3200A. This gives the results shown in 
Figure 1. 

In the temperature range we are considering, photoionization of the ground state 
of hydrogen usually dominates the total opacity chinu in the Lyman continuum. With 
T^i^u > I'Lyman) ^ 1 at the photosphere corresponding to lower frequencies, we here 
assume an LTE ground state population. We have used — uq for case I, and 
rig ~ (uq/^iY^'^ for case N, matched at 9eq. We choose uq as a typical density (in the range 
no = uh = 10^ cm~^ — 10^^ cm~^). Since this density only enters in the Lyman continuum 
it affects our analysis at longer wavelengths only through the Rosseland opacity, and as will 
be shown later this has only a very small effect on the dilution factor. Since there is little 
flux in the Lyman continuum in the temperature range we are considering, there are only 
small changes 5%) in Xn/^eO'es as we vary no through the range indicated above. 

As a first order approximation to the line scattering opacity due to (mostly iron peak) 
heavy elements, we utihze a fit to the results of WPV [lines 3 and 4 of equation (pO|)]. Since 
the strongest dependence of the scattering opacity is on temperature, the fit for the Balmer 
and Paschen continua was produced at the fiducial values of uh = 10^^ cm~^, r/v = 10 days, 
and heavy element mass fraction Z = 0.02. Varying the coefficients of the fit by a factor of 
two produced values of Xn/Xes that were at most 15% greater (doubling the coefficients) 
or at most 10% smaller (halving the coefficients). From WPV we see that changing the 
coefficients by this factor would also correspond to varying uh from 10^° — 10^^ cm~^, 
r/v = t — to from 4-30 days, or Z by a factor of 10. The largest effect these variations 
will have on our later analysis is directly through the quantity Xu/Xes at the frequencies of 
interest. 

In order to assess the importance of the absorptive part of the continuum opacity due 
to the lines, we need to estimate the thermalization parameters and the optical depths 
Tn of the lines. An estimate of the thermalization parameter (appearing in the source 
function 5"^ = (1 — e„)(Jj,)„ + e^-B^ for line n in a two-level atom without continuum) is 
provided by equation (5.3) of WPV, giving ^ 2 x 10"^ for T 5000 K, rie <^ 10^^ cm~^, 
and A ^ 6000 A (where the lines are important). The ratio of effective absorptive to 
total effective continuum opacity in the lines is then [using equations (3.10) and (3.16) 
of WPV] Klines/ Xiines ~ (en(l + Tn)) if (e„r„) ^ 1. This is not larger than the ratio of 
the hydrogen bound-free to total opacity, using the results of WPV which indicate that 
(r„) ~ 1. Therefore we shall neglect the contribution of lines (as well as photoionization of 
heavy elements) to the absorptive opacity. 

The assumption of low-density is implicitly introduced when we neglect the free-free and 
bound-free hydrogen opacities (only in our calculation of Xv and xb) for frequencies below 
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the Lyman limit. This neglect makes Xr niost about 12% too small for Uh = 10^^ cm~^. 
This implies a corresponding upper limit on uh in order that our assumption of a scattering 
dominated system be valid. As one progresses to longer wavelengths, the ratio of absorption 
to scattering increases; so our analysis becomes less accurate in the infrared. 



3.1. Thermal Structure of the Atmosphere 



Following the approach used by HLW, we next need to determine relations among the 
variables T, r and r for our two ionization regimes (I,N). A formalism for the general case 
may be developed in the following way. The derivative of equation ([T6|) with respect to r 
gives 

dz^ 16 fe,V dlnd 



Q 



(21) 



dr 3 \9 J dr 

We also have the definition (also only valid at large optical depths) 

dr^/dr = -A{rjrfn^a^s , (22) 

where A = Xn/Xes ~ A{9) has a weak dependence on the other properties of the 
photosphere: density, velocity, and heavy element abundance. In what follows, we 
sometimes take A{6) (seen in Figure 1) to be a slowly varying function, so we may neglect 
its 6 dependence in integrations involving more rapidly varying functions over restricted 
ranges of temperature. 



If we find expressions for ne{6,r) we may equate equations ( pi]) and (^2]) and obtain 
relations for 6{r/r^) in both the ionized and neutral cases. Implicit in the following is the 
assumption that near (and below) the photosphere 1 + di = rii/nl ^ 1, which we find to be 
verified by complete atmosphere calculations. 

For the ionized case rin = rie, giving 

ne = n,{rjrr. (23) 
We then obtain the temperature profile 



4 



a+l 



with Gi = y4A(^)n^,(TesAri and Ari = r^jia + 1). One expects that Gi ~ r^, so it is 
reasonable to neglect Cj ~ 1 when scattering dominance produces a photospheric optical 
depth Tp ^ 1. 
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For the neutral case uh = n\ gives ng = (n///$i)^/^, so we have 

n, = 2.94 X (r * /r)°/' ^"^/^ exp(-7^) , (25) 

where 27 = (/ii/i//c)/(5040 K) = 31.31. For this case we are not able to obtain an exact 
analytical expression. However, since we are working in the regime where 9 = 9^ = 1^ the 
integral is dominated by the exponential term and we obtain the approximate form 



(j^ 9-'''eMlO) = Cn-Gr^ 



* ^ , (26) 



r 



where Gn = 8.63 x 10^'^A{9)nJ^aesAr2 and Ara = 2r,/(a + 2). 

Since we expect the function 9{r) for the ionized and neutral regimes to match at 
rie/nH = 0.5, we determine the location 9m of the match by solving for the intersection of 
equation (p^), with Ci = and nn^i = 2: 

27^- + 1^^™ = 48-90 - Inn, + ^-(InG, - 41n^e) , (27) 

2a + 2 a + 1 

with Tm/r^ = G]/''°''^^\9m/9e)'^^''°'^^^ ■ Since \n9m = ln6'e = 0, we may approximate equation 

m as 



9m = 3.194 X 10 



-2 



48.90 - Inn, + InGi 

a + 1 



{2i 



Equation (|28|) is used only to indicate the major dependencies of 9^., and is never used in 
our numerical work. From equation (^) we find that 9^, is essentially a universal function 
of a particular combination of model parameters, indicated in Figure 2. 

Also matching equation (EH) at this point determines the integration constant to be 



//) \4 /n \ 2(a+2)/(a+l) 

Gn = C/^exp(7ft.) + G„Gr("-^^)/^^"^^M|^ , (29) 



2a + 4 

9, \ 



9.30 X loM nf^ {A.99A{9m)ARM)^^ 

\ 9rr,. I 



^ ^ 2a + 2 



■ (30) 



We have introduced the dimensionless parameters rin = n,/10^^ cm~^ and Ai?i4 = Ari/10^^ 
cm = [r*/)a + 1)]/10^'^ cm. A more aesthetically pleasing location of the match is provided 
by requiring that both the temperatures and their derivatives match. This approach yields 
a value of fm = {ne/nH)match = 0.46 — 0.51 for a wide range of models. 

In Figure 3 we plot the (inverse) temperature structure of the atmosphere with 
the choice AR14 = 1.0, riu = 1, and a = 8 (for three different values of the effective 
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temperature). As shown in this figure, these models also have an asymptotic temperature 
^max given by equations ( ]26| ) and (^9]) as 

In addition, we obtain a minimum value of the ionization fraction {ne/nH)min ~ 0.2 for 
these atmospheres with a wide range of model parameters. This validates our neglect of 
the contribution of elements other than hydrogen to the electron density. However, it is 
important to remind the reader that these temperature structures are only valid for r,^ 1. 
At small optical depths the temperature and ionization are governed by the radiation field 
(formed at large optical depths), which tends to drive the level populations far from their 
LTE values. 



1 + 



+ 1 

2V2 \a + 2 



(31) 



3.2. Location of the Photosphere 



Following the approach of HLW, we now operate on our master equation (|T^) with 



V ^du, obtaining 



^ qv 



A 



ff 



i-1 



D, + J2A[[Di-ElJ2Mr,'D,/E', 



1=1 



k=l 



l=j \ k=l J 



Ues fc 



(32) 



where 



/■CXD 

= / (J.-i?.)(l-e-'^^/'=^)-V-4cii 



El 



{c^/2h) 



and A'l = Ai/Aff. Next, we invert the above matrix equation to obtain the first-order 
quantities and substitute the expression into the equivalent quantities in equation ([T3|). 



We next approximate the resulting transfer equation as the photosphere is approached 
from below by replacing all occurrences of J,y with B^, but keeping terms in J,, — By as the 
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first order correction. This yields (for z/„ < z/ < 



3/1 



l=n j=l ^3 



witfi 



1 



dr, 



(33) 



(34) 



qu 



and witli tlie integral in (dimensionless) E'l now over instead of J^. Since all the Gaunt 
factors are unity, 

5.20 X 10^6^^^cmh~^ , 



A 



ff 



A'l = 62.51-^6 exp{31.3ie/ 1') . 



(35) 



The structure of the matrices Qij and Mij depends upon the approximations employed. 
Naively, in the limit of very low densities collisions may be neglected, giving Mij = 6ij, as 
in HLW. However, since the Lyman continuum is usually in radiative detailed balance due 
to its large absorptive opacity, we keep the coUisional coupling between the n = 1 and 2 
levels of hydrogen. This modifies M^^ by the inclusion of the term u = C'y^l (1 + C'y^ in 
two elements, M^x = \ — u and M{^^ = m, with the rest of the matrix being unchanged. 
Because of the relative lack of photons in the Lyman continuum at the temperatures of 
interest, u is not necessarily small, even at our low densities. We have explored the range 
< n < 1, and the results we report for the optical and IR do not change as we vary u. 
The largest change occurs at the highest temperatures (T ^ 12,000 K), but even then the 
changes are relatively small and are restricted to the Lyman continuum. 

We now investigate the properties of the photosphere in terms of the major model 
parameters: 6'e, A_Ri4, rtl^, and a. In keeping with our first-order analysis, we may use 
the zeroth order expression (|16]) for the relation between optical depth and temperature in 
simplifying the expression for d^By/dr'^^ . We obtain 

d^Bjdrl = i9/lQ)ixR/x.fiO/e,)'BJihu/kT) , (36) 

where 

/(-) = '""'r/^^if'" . (37) 



16(1 



The first order deviation of the average intensity from blackbody is then given by 



\Ji, — Bu 

b7~ 



2.778 X 10 



Uf, cm'^ 



10 



4 5)9/2 



(3^ 
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in the interval z/„ < i/ < z/„_i, with 
SiiO,u) = xr/Xv , 

= XR/Xes, (39) 

/-co 

Qn = E^^E(^^7V^;)EQ7fc / 5i(^,^)/(a:)(e^-l)-Vda;, (40) 

i=n 3=1 k=l 

and ?/ = hv/kT. In order to further utihze these equations, we must determine the function 
?T,~^(r*/r)^, which is straightforward in the two hmits of ionization. 

For the ionized case (I) we may write equation (^) in the following form after using 
equations (p3D and (|2^ ) to obtain rieiO) and r{0) : 



\J,-B,\/B, = Fi{ARi^,nn,a)F2{ee,a)F^{e,v,a) , (41) 



where 



(a-4)/(a+l) 
14 



Fi = (4.99nn)-^/("+^)Ai? 

In the case where Mf^-^ = 6ij (no collisions), a ^ oo (sharp atmosphere), and Xu = '"-eO'es 
(no line scattering), this case corresponds to the treatment of HLW. 

We estimate the location of the photosphere at each frequency as the depth at which 
the mean intensity differs from the Planck function by of order the Planck function, 
\J,^ — B^\ ^ Bi,. As in HLW, we fit a smooth function through the temperature region where 
Ji, — Bu changes sign. This procedure should reflect the smooth dependence of the location 
of the photosphere on our model parameters. Numerical calculation of Si, S2, and S3 
allows us to determine the location of the photosphere, represented by 6p (for fixed z/ and 
a), as the particular function of model parameters 9e, Ai?i4, and rin indicated in the 
above functions Fi and F2 and shown in Figure 4a. 



The neutral case (N) is much more difficult. Inspection of equations ( PBj ) and 
shows us that we cannot separate the dependencies as completely as in equation (|4T|) since 
model parameters occur in both C„ and Nevertheless, using equations (|^) and (^), 
equation (|38|) now assumes the form 



B, V 'a + 2 



a-S -10 



- 10.85 (aRu^) f3.63 x lO^^i) 9; 
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X sl^+' e ^ e^' (Cn -ele^ e^') ' 53 . (43) 

Figure 4b is a plot of the numerically calculated 6p [from setting equation (^31) equal to 
unity] for a range of model parameters, with the constraint that 6'^ < < Omax- If Tp is 
only slightly less than T^, then Op should predominantly be a function of the combination 
X = ^(i0"-8)/{a+i)^^(8-")/(2"+2)^9/(2a+2)^ ^^^.^^ ^^^^ ^^^^^^ Special case of a = 8, 

there is some spread in the relation. However, for a = 6 we obtain a narrow spread in the 
relation 9p vs. X, corresponding to 9p/9m < l-l- The spread is much greater for a = 10. 



3.3. Calculation of the Dilution Factor 



Now we have all the ingredients necessary to calculate the dilution factor ( from 
equation (|1^) in terms of model parameters. In the previous section we indicated how the 
function 9p{9e, ARi4^,nii,a; v) is obtained. Thus for any desired set of frequencies we may 
replace 9e through this function in equation ([T^), which yields C = C(^p! ""-ii) -^-^14! 



Thus far, the model parameters have referred to a fiducial radius r^, and thus are not 
observed quantities. We now calculate the dilution factor in terms of parameters evaluated 
at the photosphere, which are (potentially) observable. In addition to 9p{v) we employ 
the photospheric radius Vpiy) = Vp(t — to), where the velocity Vpi^u) at the photosphere is 
obtained from analysis of an appropriate line profile (most reliably from the sharp minimum 
of a weak line). The photospheric density is then n/^(rp) = np{v). 



We now choose our fiducial radius to be the photospheric radius, so that 



V) = rpl4 X 10^^ cm and = np{v) = Upn x 10^^ cm We also define 



the photospheric scale height Ai?pi4 = rpi4/{a + 1). Setting equation (|38D equal to unity, 
we then obtain relations among the photospheric quantities: 

Q.2779n;,\9:'9l/^S2i9p)Ssi9p) = 1 (44) 

for the ionized case and 

2.988 X W-\li^9:'9f/^e^''^S2{9p)S3i9p, u) = 1 (45) 

for the neutral case. We use these two equations to eliminate n^n from equations (p^ ) and 
(p6|) at r = = r^. 

We also use equation (^) to eliminate 9^. in favor of yielding expressions for the 
dilution factor for both cases in terms of photospheric quantities. For the ionized case, we 
obtain 

C = = QM7AR-pU'el"'Sl'p'S^p"'S^p"' , (46) 
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and for the neutral case we obtain, for ( = 

'a + l\ 1.96 



C '^^RpU^l^'^S^pS2pS^p + 10.85 (^^ ^ 2^ S2p - S2m 



10.85 



a + 2/ e. 



a + 2 



1.498 X 10^C52"^^\-'A^r4 /'^i'/'e-^^^Sf S,-/^3-/] , (47) 

where the subscripts p and m indicate that the function is evaluated at 6p and 6m, 
respectively. These equations represent the most important results of this paper. The 
dilution factors for the mostly ionized and mostly neutral cases are presented in Figures 5a 
and 5b. 

For the ionized case we have a universal dependence on photospheric scale height, with 
a unique dependence on photospheric temperature for each wavelength, as shown in Figure 
5a. The decrease in the dilution factor (at fixed scale height) with decreasing temperature 
and wavelength is due to the increase in UV opacity as the temperature drops, although 
the (Rosseland) mean opacity is barely changing. 

For the neutral case there is no such simple scaling. In order to obtain the dilution 
factor from equation ( ^7] ) we employ the same steps as above to eliminate n^, Ari, and 9^ 
in equation (^7]) for 9m in favor of (, ARpi^, and 9p. We then iteratively solve equation 



(W^), which is easily done because 9m is a weak function of the parameters, as indicated in 
Figure 2. In addition, we only accept solutions for which 9m < 9p < 9max- For a > 4 it 
is found that there are either two or no solutions for The smaller root always occurs 
where 9p/9max ^ 0.99, which corresponds to the outer crossing of the rie/uH = 0.5 and 
temperature curves, as seen in Figure 3. We see that the physical location of this root is 
in the outer region of the atmosphere that is not accurately modelled. In addition, this 
smaller root does not match Ci as the photosphere becomes ionized. For these reasons we 
always choose the larger of the two values of Cn- 

Figure 5b shows how Cn varies with photospheric temperature, scale height, and 
wavelength; and its match to Q. We also find that for 9p only slightly greater than 9m 
we obtain a dependence Cn oc Ai?"// as expected when the first term in equation (^) 
is negligible. For fixed scale height ARpn = Vpi^/^a + 1), the behavior of the dilution 
factor shown in Figure 5b is relatively insensitive to the value of a (at least for the range 
6 < « < 10 investigated). 

A particular frequency dependence that we have investigated is that across the 
Balmer ionization threshold at A = 3646 A, motivated by the dependence on AR found by 
Hershkowitz, hinder, & Wagoner (1986a). We find that the absolute value of the fractional 
jump in the photospheric temperature (corresponding to the results in Figures 4a and 4b) 
is less than 0.03 for our ranges of the photospheric parameters. The fractional jump in the 
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dilution factor (corresponding to the results in Figures 5a and 5b) is likewise found to be 
less than 0.04. Therefore it appears that the line scattering (and atmospheric extension) 
has reduced the jump from the values found with only electron scattering in a sharp 
photosphere by Hershkowitz, Linder, & Wagoner (1986a). 
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4. Discussion 



The determination of the distance of a supernova via EPM would proceed as follows 
within the above formulation, which assumes nothing (except spherical symmetry) about 
the nature of its progenitor. (The other assumptions we have listed in section 2 can be 
checked after one obtains the photospheric conditions from our analysis.) From each 
observed spectrum, the photospheric temperature Tp(z/) is first estimated by fitting Planck 
functions to the continuum in the neighborhood of various frequencies. As has been 
indicated above, the radius rp(z/) of the photosphere is obtained from (the weaker) line 
profiles. The remaining parameter to be determined is a, the total hydrogen density 
power-law index. 

A comparison of detailed model atmospheres with spectra of SN1987A obtained during 
days 2-10 after explosion led Eastman & Kirshner (1989) to conclude that 7 < a < 11. 
It was found that the UV continuum as well as the line shapes were sensitive to this 
parameter, although not greatly so in the range indicated. However, Branch (1980) has 
shown that the effects of optical depth (i.e., heavy element abundance) and density profile 
on line shapes may be difficult to separate. In addition, nonLTE effects (which he did not 
include) will be important for some lines. 

Since it may be difficult to determine a accurately from lines, let us consider whether 
it might be obtained in another way from observations of the continuum. If one follows a 
fiducial volume element which always contains the same nuclei, its radius and density soon 
obey oc t — to and oc (t — to)~^- K then follows that the corresponding photospheric 
quantities are related by n^iy.t) oc (t — t^'^°'~'^^r~'^{y^t). From this relation, one then sees 
that 



(x{t) 



9 In n„ 

O ~r 



ain(t - to 



1 



d\n{t-to 



(4J 



In principle, one could obtain the function np{9p,rp,a) by using equation (|T9|) to write 
equations ( ^^ and (|1^) in the form 

Up = 2.78 X lO'\t9;'/''S^^^S2pSsp cm'^ (49) 

for the ionized case and 

Up = 0.893C0;''^'e'^''S^p'SlpSl cm'' (50) 
for the neutral case. Our results for the dilution factor C{9p, Vp, a; v) would then be inserted 

j-2/3^_7/6r,-l/3 ^^1/3 
t.pl4 Op D2p O^p 



into these relations. For instance, one would then obtain riu = 0.224Ai? ]^^^6' "^^^82}^^ Sl^^ 
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for the ionized case. We also note that SKE quoted the same relation between ( and Hp 
as seen in equation (^91), for fixed 9p. Of course, since Up depends upon a, the solution of 
equation (^) would require iteration, using values of 9p and Vp obtained at various epochs. 
One could assume that a{t) was a slowly varying function. However, a test of the practical 
viability of this method is beyond the scope of this paper. 

Once the dilution factor C,{9p,rp,a]v) has been determined, its relation to the 
luminosity [equation (|18|), with = r^] could be applied in two steps. First, because of 
the frequency dependence of the dilution factor, the previous estimate of the photospheric 
temperature should be improved by iterating the fit of equation ( [T8| ) to the observed 
continuum. Second, once the fundamental parameters Tp(z/), rpi^u), and a have been 
determined, equation can be employed to obtain the luminosity distance. Of course, if 
the supernova is in the Hubble fiow, the redshift Z (of the supernova or parent galaxy) will 
then produce a value of the Hubble constant [and for redshifts Z > 0.4, the deceleration 
parameter (Wagoner, 1977)]. 

Another important question is the sensitivity of the dilution factor to the coefficients 
of the line-scattering opacity. We have both increased and decreased the coefficients by 
a factor of two. The dilution factor is most affected in the Balmer continuum (where 
the effects of the line opacity are strongest). At the highest temperatures, the dilution 
factors in the Paschen continuum only vary by a few percent. However, as the scattering 
increase at lower temperatures, we see much larger affects. Equations (^6]) and (|39D show 
that Ci oc 5*1^ , so doubling the coefficients leads to variations of about 25% when the 
line scattering is greater than the electron scattering. However, doubling the coefficients 
corresponds to increasing the heavy element abundance by about an order of magnitude, as 
shown by WPV. In principle, this abundance can be determined by the behavior of the UV 
portion of the spectrum. 

Like the heavy element abundance, the extinction of the spectrum by our galaxy and 
the supernova parent galaxy must be determined before a reliable luminosity and distance is 
obtained. If the dependence of the extinction on wavelength is universal, then its magnitude 
can be obtained by including this dependence in the fit of the spectrum. However, this 
procedure only becomes reliable and sensitive if the observations extend from the IR to at 
least the near UV. 

A major challenge that faces us is to reconcile our dilution factor with that obtained 
by SKE. There is qualitative agreement in the ionized regime, although the temperature 
dependence is somewhat different. However, in the recombination era our dilution factor 
decreases as the temperature decreases, whereas the SKE dilution factor [as well as that 
recently obtained by Baron, et al. (1994)] increases. We can understand our result 
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based on the fact that as the ionization fraction decreases, the opacity due to absorption 
(proportional to n^) decreases faster than the electron and line scattering opacity (roughly 
proportional to ng), producing a more dilute radiation field. Another point to make in 
comparing our work to that of SKE is that the color temperature they determined may 
not correspond to our photospheric temperature. These temperatures should be similar as 
long as there is no net flux in the lines, the frequency dependence of the dilution factor 
is negligible, and a continuum can be uniquely determined. The increase in the density 
of lines toward the UV make the determination of a photospheric temperature potentially 
more difficult in the B and U bands. 

Some other crucial questions remain, which can only be answered by detailed 
comparisons between the spectral luminosity predicted by this method and the observed 
flux of a variety of Type II supernovae at various epochs. This is the next step in our 
program. Some of these questions are: 

1) How valid are the approximations (such as scattering dominance) that we have made in 
this analysis? (If the dilution factor approaches unity, as probably occurs in the infrared 
wavelengths, the photosphere is no longer scattering dominated.) 

2) How closely does our definition of the photosphere correspond to reality? 

3) How tightly can the luminosity be determined from the ranges of the parameters 
constrained by the spectral fit? 

In spite of these uncertainties, we believe we have developed a potentially useful 
tool for determining the luminosity of Type II supernovae directly from observables in a 
model-independent manner. While the results for the ionized regime seem fairly robust, 
more work is needed to understand more fully the applicability of our results in the neutral 
regime. 
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Fig. 1. — The temperature dependence of the ratio of Rosseland mean to total opacity, 
choosing fiducial parameters Uh = 10^^ cm~^, r/v = 10 days, and Z = 0.02. The curves 
correspond to the following wavelengths, from top to bottom: 9000 A, 7000 A, 5500 A, 4400 
A, 3650 A, and 3000 A. The top curve also gives Xr/Xcs, since for A > Xpaschen-, Xiy ~ Xes- 

Fig. 2. — The (tight) dependence of the matching temperature 9m on a particular 
combination of model parameters. This plot has points for a = 6, 8, and 10. The range of 
the other parameters is 0.01 < Ai^M < 100 , 0.01 < nil < 10, 0.6615 < 6^ < 1.26. 

Fig. 3. — The (inverse) temperature profile as a function of the scaled radius for model 
atmospheres with a — 8, riu — 1 and AR^ = 1; for several values of effective temperature. 
Note that the curves match smoothly at the fitting point Ue/nn = 0.5. The ionization fraction 
remains greater than 10% for a wide range of models, temperatures. 

Fig. 4a. — The dependence of the photospheric temperature, in the ionized case, on a 
combination of model parameters for three choices of wavelength. The dashed and solid 
curves are for a = 10 and a = 6, respectively. The results for A = 3650 A and 5500 A are 
essentially identical to that for 4400 A. 

Fig. 4b. — The dependence of the photospheric temperature, in the neutral case, on a 
combination of model parameters for A = 5500A. The ranges of model parameters are the 
same as in Figure 2. There is no photospheric solution when the abscissa is less than about 
— 1.5 (representing low densities and high temperatures, the ionized case). 

Fig. 5a. — The (scaled) dilution factor for the ionized case as a function of photospheric 
temperature. The behavior at short wavelengths refiects that of XrIXv- 

Fig. 5b. — The dilution factor as a function of photospheric temperature for a = 8, at three 
wavelengths. Both the ionized regime (dashed) and the neutral regime (solid) are shown. For 
each wavelength, the upper curve is for Ai?pi4 = 1 while the lower curve is for Ai?pi4 = 10. 



